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ABSTRACT 

Carmichael C 1 3 has proposed a direct injection, 
temperature ignition <TI), r ec i pr ocat i ng engine design which 
divides the conventional four stroke cycle functions between 
a compression cylinder and an expansion cylinder, 
interconnected by a static regenerative heat exchanger. The 
prediction of heat transfer loss in the expansion cylinder is 
required to predict the performance of the new cycle. 

Woschni's correlations [93 , developed for direct 
injection, quiescent, four stroke, compression ignition (Cl) 
engines, are used to predict uncorrected expansion cylinder 
heat transfer loss. Laws of similarity allow, in principle, 
extrapolation outside the experimental data range for forced 
convection heat transfer only. General applicability of 
Woschni's correlations is established by examining the 

expansion cylinder quiescent condition, using predicted mean 
gas velocity magnitudes and fully developed, turbulent flow, 
forced convection models. Predicted radiant peak flame 
emissivity and temperature are found and compared with 
nominal peak values observed in conventional Cl engine cycles 
to establish a minimum radiation flux reduction ratio. A 
mean reduction ratio for forced convection flux is also 
found. A corrected, conservative prediction of heat transfer 
loss in the expansion cylinder is then established. 

Thesis Supervisor: A. Douglas Carmichael 

Title: Professor of Power Engineering 
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NOMENCLATURE 


A1 phanumer i c 




A 


heat transfer surface area Cft' 1 '] 


a ht 


cylinder heat transfer surface area Cft^l 


a p 


regenerator port exit cross-secti onai area 
Cf t 2 ] 


cl 


crank radius length Cft3 


B 


bore diameter Cft] 


B e 


fuel consumption rate Clbm/hr] 


Bo 


Boltzmann number; defined by equation (61) 


c -f 


friction coefficient; defined by equation (11) 


s 


specific heat at constant pressure 
CBTU/lbm— °R] 


C v • 


specific heat at constant volume 
CBTU/lbm— °R 3 


c x & c 2 


Woschni correlation constants; defined by 
equations 18a - lSe 


c 


speed of light = 2.997902 x 10 ^ cm/sec 


c b 


Boltzmann's constant = 1.38026 x' 10 - *k erg/°K 


C P 


Planck's constant = 6.623777 x 10 ~^ erg— sec 


c s 


soot formation rate coefficient 
C 1 bm/sec- f t°— atm] 


D 


diameter Cft] 


E 


"radiant emissive power CBTU/hr-f t^ or 
BTU/sec-f t 2 ] 


' E ;v 


radiant monochr omat i c emissive power 
Cergs/sec-cm -A (cm) ] 


E s 


soot combustion activation energy 
CBTU/mole (lbm) ] 


F 


radiant geometric view factor 


h 


heat transfer coefficient CBTU/hr-f t^-° R] 
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h 



D 



h 



e 



h 



1 



j 

K 



k 



f 



k 



g 



k 



p 



k 



A 



L 

L.0 

Z 

1 



1 



c 



1 



R 



1 



r 



m 




m 



N 



RF'M 



Nu 



P 



heat transfer coefficient across pipe^ 
diameter D; turbulent flow CBTU/hr-f t^— °R] 

total cylinder heat transfer coefficient; 
turbulent flow C BTU/hr-f t 2 - 0 RD 

heat transfer coefficient across length 1 of 
flat plate; turbulent flow CBTU/hr-f t - * 1 - 0 RD 

local heat transfer coefficient at distance x 
from leading edge of flat plate; turbulent 
flow CBTU/hr-f t^- 0 R] 

Colburn factor; defined by equation (9) 

thermal conductivity CBTU/hr-f t-° R3 

apparent flame absorption coefficient per 
unit length Cft - *! 

apparent gas absorption coefficient per unit 
length Cft -1 ! 

apparent flame absorption coefficient per 
unit length-unit pressure Cf t-atml - * 

monochromat i c absorption coefficient per unit 
length Cft-A3 - '*‘ 

piston stoke CftD 

cylinder gas volume height CftH 

flame luminosity 

flat plate length Cftl 

Re char acter i st i c length Cftl 

effective radiant mean beam length CftD 

connecting rod length Cftl 

charge mass Clbml 

charge mass at ignition start ClbmJ 
charge mass flow rate Clbm/secl 
cylinder cycle frequency CRF'Ml 
Nusselt number; defined by equation (3) 
cylinder pressure Catml 
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p 



c 



p 



o 



Pr 



P 



u 



Q 



Q 

q 



R 



c 



R 



o 



R 



Re 

Re 



D 



Re. 



r _ 
*_ 




.3 

SFC 

St 



s 



T 

Tg 

T 3c 



T 



o 



cylinder pressure at ignition start Catml 

cylinder motoring pressure Catml 

Prandtl number; defined by equation (5) 

unburnt -fuel partial pressure in the soot 
-formation zone Catml 

heat transfer CBTU3 

heat transfer rate CBTU/sec] 

heat transfer flux CBTU/sec-f t^H 

monochromatic radiant heat transfer flux 
CBTU/sec-f t 2 -x 1 

crank ratio = l r /a 

universal gas constant CBTU/mol e ( 1 bm) RJ 

specific gas constant = 53.90 lbf-f t/1 bm-°R 

Reynolds number; defined by equation (4) 

Reynolds number for pipe flow; 1 = D 

local Reynolds number for flat plate at 
distance x; 1 = x 

cylinder compression ratio 

piston speed Cft/sec3 

soot concentration Clbm/ft"'] 

specific fuel consumption Clbm fuel /HP— hr D 

Stanton number; defined by equation (10) 

straight line distance between cylinder crank 
shaft and piston pin Cftl 

temperature C°R3 

cylinder bulk gas temperature C°R3 

cylinder bulk gas temperature at ignition 
start C°R3 



cylinder motoring bulk gas temperature C°R3 



T u 


local temperature of the soot -formation zone 

[or: 


T w 


heat transfer surface wall temperature C°R3 


t 


unit time CsecU 


V 


cylinder gas volume C f t □ 


Vc 


“jr 

cylinder clearance volume at piston TDC C f t 3 


V c 


cylinder volume at ignition start (ft"’) 


V u 


•*7* 

volume of the soot formation zone Cft'"’3 


U 


bulk fluid (gas) velocity Cft/secH 


Ui 


impulse charge velocity exiting regenerator 
port Cft/secH 


Uex p 


bulk gas velocity at the cylinder walls in 
the expansion stroke [ft/sec3 


Uex h 


bulk gas velocity at the cylinder walls in 
the exhaust stroke II ft /sec 3 


U> 

r 


angular gas velocity at radius r Cft/sec3 


X 


radiant path length Cft3 



Symbol s 



§ 


charge equivalence ratio 


$ 

*u 


unburnt fuel equivalence ratio in the soot 
formation zone 


© 


crank angle C° degrees3 


S 


forced flow boundary layer Cft3 


A 


electromagnetic wavelength Ccm or /un3 


A 


dynamic viscosity II 1 bm/f t-sec 3 


P 


density Clbm/ft"‘3 


V k 


kinematic viscosity (momentum diffusivity) = 
}i./p Cft j:: ’/sec3 


« 

r 


thermal diffusivity = K/C^p Cft" i /sec3 
ratio of specific heats = Cp/C v 



T 



w 



£ 



£ 



a 



£ 



f 



£ 



g 



£ 



o 



£ w 
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*w 



A 

A heat 



tl 

y 



apparent shear stress acting on the -fluid at 
the wall Clbf/ft^D 

emissivity coefficient; grey body 

apparent -flame emissivity coefficient; grey 
body 

■flame emissivity coefficient; grey body 

gas emissivity coefficient; grey body 

emissivity coefficient of infinitely thick 
radiating flame; grey body 

wall emissivity coefficient; grey body 

transmi ttance 

absorptance 

monochromatic absorptance 

ref 1 ectance 

wall reflectance 

—9 

Stef an-Bol tzmann constant = 1.713 x 10 
BTU/hr-f t 2 -°R 4 

incremental change or difference 

specific heat release per cycle 
C BTU/1 bm— char ge ] 

heat transfer loss/heat release per cycle 

nominal Qp/Qj per cycle in conventional 
four stroke, high r , direct injection Cl 
engi nes 



Subscr i pts 



cp 


combustion products 


FC 


-forced convection heat transfer 


NCE 


New Cycle Engine 


R 


radiation heat transfer 


T 


total heat transfer 


W 


Woschni test engine 


e 


crank angle degrees 


Superscr i pts 





mean value 



* 


corrected value 


Abbreviations 

BDC 


bottom-dead -center 


TDC 


top— dead -center 


BTDC 


before TDC 


C/H 


carbon/hydrogen mass ratio 


Cl 


compression ignition 


SI 


spark ignition 


T I 


temperature ignition 



12 





LIST OF FIGURES 


Figure 


Title Page 


1-1 


New Cycle Engine vertical cross- 21 

section C13 


1-2 


New Cycle Engine expansion cylinder 22 

and piston linkage dimensions C 1 3 


1-3 


New Cycle Engine firing cycle P-V 23 

diagram; natural aspiration Cl 3 


1-4 


New Cycle Engine firing cycle Tg-V 24 

diagram; natural aspiration Cl 3 


2-1 


Woschni 's test engine firing cycle; 41 

natural aspiration C93 


CM 

1 

CM 


Expansion cylinder mass charge, m, 46 

versus expansion stroke crank angle 


2-3 


Expansion cylinder charge Ui versus 47 

expansion stroke crank angle 


3-1 


Expansion cylinder firing cycle h 0 62 

and Q versus crank angl e; Woschni ' s 
correl ati on 


4-1 


Qp/Qy for two 4 stroke, direct 73 

injection Cl engines under varying 
operating condition as observed by 
Qguri et al C83 


4-2 


Soot production fuel jet mixing 75 

model of Khan et al C183 


4-3 


Comparison of in-cylinder soot 78 

concentrati on , .3, versus crank 

angle as observed by various 

investigators C103 


4-4 


Variation of kp versus angular 83 

displacement for £ a calculation; 

Sitkei et al model C143 


4-5 

4-6 


Measured T^ in Cl engines versus 85 

crank angle as observed by various 
investigators C 103 

Flame luminosity, <£, as a function 87 

of fuel C/H mass ratio in low 
combustion pressure power units C193 




13 



4-6 



87 



List o-f Figures (Cont) 



F i gure 


Title Page 


4-7 


Variation o-f Tg/T^ in the function 89 

of the Bo number C143 


A— 1 


Generation of expansion cylinder 108 

firing cycle motoring P Q 


B-l 


Standard reci procat i ng cylinder and 114 

piston linkage geometry 


C-l 


Ratio of specific heats, Jf, of lean 117 

hydrocarbon-ai r combustion products 
as a function of P, Tg, and §; fuel 

C n H 2n C73 


C-2 


Cp of lean hydr ocar bon-ai r 118 

combustion products as a function 

of P, Tg, and $; fuel CH~ C73 

1 1 «•£.! 1 



14 



LIST OF TABLES 



Table 



2-1 



3-1 



A— 1 



C-l 



Title Page 

Calculation of order of magnitude 40 

reduction in combustion forced 

convection q from Woschni's test 

engine to New Cycle Engine 

expansion cylinder operation at 5° 

crank angle increments 

Calculation of charge Ui leaving 45 

the regenerator port into the New 
Cycle Engine expansion cylinder at 
5° crank angle increments 

Comparison of h g calculated by 52 

Woschni's correlation and 
theoretical models during charge 
induction into the New Cycle Engine 
expansion cylinder 

Calculation of New Cycle Engine 58 

firing cycle Q using Woschni's 
correlation at 5° crank angle 
i ncrements 

New Cycle Engine expansion cylinder 106 

computer-generated thermodynami c 
firing cycle model as a function of 
crank angle; 5° increments Cl] 

Coefficients for calculation of Cp 115 

and ji. of principle gases as a 
function of Tg <°R) C6] 



15 



CHAPTER ONE 



INTRODUCTION 



1 . SCQF'E 

Conventional compression ignition (Cl) reci procat 1 ng 
engines have -found significant application in commercial, 
industrial, and in land and marine transportation systems. 
Recent emphasis on -fuel costs and avai 1 abi 1 i ty , in 
conjunction with improved materials engineering, has prompted 
consideration of design alternatives which offer increased 
thermal efficiency and reduction in operating costs. A New 
Cycle Engine design has been proposed by Carmichael C 1 3 in 
which the functions of a conventional four stroke 
reciprocating cylinder are divided between a compression 
cylinder for charge induction and compression, an expansion 
cylinder for charge combustion-expansion and exhaust, 
i nterconnec ted by a static regenerative heat exchanger. 

The attract i veness of such an engine is possible 
because improved fabrication techniques presently allow 
construction of strong, thermal resistant, extremely thin- 
walled ceramic matrices possessing high heat transfer 
properties. A cyclic interchange of compressed charge into 
and expired exhaust gas out of the expansion cylinder via the 
regenerator establishes an effective medium for heat transfer 
and allows use of exhaust heat otherwise lost to the system. 
The temperature of charge air heated by adiabatic compression 
is boosted from the regenerator heat exchange, a direct 
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consequence of which is significant reduction in cylinder 
coiiipressi on ratio, r c , required to achieve optimum 
efficiency. A computer-generated ther modynami c model of this 
engine in preliminary design indicates that improved thermal 
efficiency may be achieved. An important element in the 
cycle analysis is accurate estimation of heat transferred 
from the gas medium in the expansion cylinder to its 
enclosing surface area, and ultimately to the coolant. The 
magnitude of this heat transfer reduces both heat 
availability to the regenerator and the thermal efficiency. 



The purpose of this thesis is as follows: 



1. To evaluate the appl i cabi 1 i ty of existing 
heat transfer correlations developed to 
predict total mean instantaneous heat 
transfer coefficient, h e , in four stroke 
direct injection, quiescent, Cl engines 
for prediction of total heat transfer in 
the New Cycle Engine expansion cylinder. 

2. To evaluate and compare the magnitude of 
the reduction ratio in both the forced 
convective and the radiant heat transfer 



from the high peak combustion 
conditions inherent in high 
engines to low peak combustion 
conditions inherent in the low 
Cycle Engine design. 



pressure 
r c , Cl 
pressure 



r c , New 
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To use the reduction ratios to correct 



existing heat transfer cor rel at 1 ons , if 
applicable, and establish a correction 
allowing better prediction of heat 
transfer loss in the New Cycle Engine 
expansion cylinder. 

2. NEW CYCLE ENGINE DESIGN BACKGROUND 

A vertical cross-sect i on of the New Cycle Engine design 
is provided in Figure 1-1. For simplicity, each cylinder is 
considered of simple geometry, with strai ght • forward walls 
and flat cylinder head/piston crown con-figuration; bore 
diameter, B, equals stroke length, L. Gas velocities are 
boosted by piston speed, induction suction, charge impu'lse 
exiting the regenerator port , % and combustion turbulence only; 
quiescent flow should predominate with no deliberate design 
feature incorporated to accelerate any rotational velocity 
incurred by off-set regenerator port impulse. The cylinder 
crank shafts will be mechanically linked to enable the 
compression piston to cycle at twice the speed of the 
expansion piston, alternating between an induction- 
compression and a free motoring stroke. Proper crank angle 
phase shift will exist to sequence the compression piston 
motion at the start of expansion piston down-stroke to 
optimize charge transfer. 

Correct value timing is imperative and is sequenced 
through one charge cycle as follows: 



IS 



1 . 

o 

» 


Compression cylinder induction stroke: 
induction inlet valve 1 opens; 
regenerator transfer valve 2 is 
cl osed . 

Compression cylinder compression 
stroke: valve 1 closes; valve 2 
remains closed. 

Expansion cylinder combustion- 
expansion stroke: regenerator exhaust 
value 3 closes; valve 2 opens inducing 
compressed charge impulse through the 
regenerator as the expansion piston 


• 


starts the down-stroke; valve 2 then 

closes. 


4. 


Expansion cylinder exhaust stroke: 
valve 3 opens, regenerating heat 
exchange capacity during exhaust gas 
expiration. Valve 2 remains closed; 
valve 1 is cycling to support the next 
charge cycle as expansion piston 
reaches top-dead-center (TDC) . 


A 


vertical cross-sect i on schematic of the expansion 


cy 1 i nder 


and piston linkage showing principle dimensions is 


provided 


in Figure 1-2. A derivation of geometric 



relationships required to calculate instantaneous expansion 
cylinder heat transfer surface area, a hT’ volume, V, and 
piston speed, Sp , is enclosed as Appendix B. Instantaneous 



expansion cylinder pressure, P , mean gas temperature, Tg, and 
volume, V, at 5° crank angle increments -from the computer- 
generated thermodynamic model is inclosed as Appendix A, 
Table A-l Cl associated P-V and Tg-V diagrams for the 
entire cycle are illustrated in Figures 1-3 and 1-4 
respectively. Other principle design and operating 
parameters for the expansion cylinder required to support the 
heat transfer calculations are listed in Appendix A. 
Appendices A and B provide the data point basis for the heat 
transfer calculations. A consideration of piston ring land 
area, cylinder sleeve clearances, and clearance volume, Vc , 
at TDC is neglected. All the area and volumetric 
calculations are based on simple cylinder geometry. Within 
the accuracy of reci pr ocat i ng engine heat transfer 
correlations in general , errors associated with these 
assumptions are considered to be insignificant. 



VAI VES 

1 INDUCTION inlet 

2 regenerator transfer 

3 REGENERATOR EXHAUST 



COMPRESSION 

CYLINDER 




LINKAGE 



FIGURE 1-1 



New Cycle Engine vertical cross-section 



D1 
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BORE DIAMETER 




FIGUSi 1-2: lieu Cycle Engine expansion cylinder and oiston 
linkage dimensions Cl] 
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FIGURE 1-3: New Cycle Engine firing cycle P-V diagram; 
natural aspiration CO 
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FIGURE 1-4: New Cycle Engine firing cycle Tg-V diagram 
natural aspiration Ell 
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CHAPTER TWO 



FORCED CONVECTION HEAT TRANSFER IN THE 
NEW CYCLE ENGINE EXPANSION CYLINDER 

1 . THEORETICAL BACKGROUND 

When a temperature difference exists between a solid 
surface and a moving fluid, thermal energy is transferred by 
the joint process of conduction and convection. Shear 
stresses acting on the fluid will slow the fluid immediately 
adjacent to the surface, forming a thin boundary layer of 
thickness &. The heat transferred between the bulk fluid and 
the body is transported by conduction in the boundary layer 
as a function of fluid movement source (natural or forced) , 
flow character (laminar or turbulent), physical properties of 
the fluid, and body geometry and roughness. The heat 
transfer rate, Q, by conduction in any direction x at any 
instant is determined using Fourier's Law by a heat 
conduction coefficient, K, and the temperature gradient: 

Q = -KAVT (1) 

ox 

In laminar flow, fluid particles flowing parallel to 
the body surface conduct heat in a normal direction. At a 
certain transition velocity, turbulence is created where 
particles have a velocity component normal to the bulk stream 
direction; the intense mixing of individual layers promotes 
rapid transpor tat i on of heat to/from the surface. Even in 
turbulent flow, however, a thin boundary layer exists due to 
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fluid viscosity across which laminar flow and heat transfer 
characteristics are maintained. In turbulent systems, the 
heat transferred between a body surface and moving fluid 
remains limited by boundary layer thickness, &, and fluid 
thermal conductivity, K. 

An ideal linear velocity and temperature gradient 
across the boundary layer between a gas medium and a body 
surface allows approximation of the temperature gradient as 
(Tg-Tw)/6; substitution into equation (1) yields the 
Newtonian expression for steady flow conditions: 

Q = -KAoT 
~o x 

= -KAC Tq-Tw H 

& 

= -K(Tg-Tw) 

o 

where letting h ss -K : 

& 

Q = hA(Tg-Tw) (2) 

For real boundary layers of varying velocity and temperature 
gradients, h may be determined in some cases by integration 
of the differential equations describing the heat transfer 
process. For more general and practical app 1 i cat i ons , 

however, laws of similarity may be used to predict the value 
of h. According to the laws of similarity, it is sufficient 
to determine the character i st i c relations between independent 
dimensionless groups obtained, recognizing that any such 
relations will apply to all similar heat transfer processes. 



From the Buckingham pi theorem which allows reduction of 
initial variables to minimum dimensionless groups, the 
following dimensionless numbers are used to describe forced 
convection heat transfer in a fluid medium to a solid body: 



1. Nusselt #, Nu = h 1 c / K (3) 

2. Reynolds #, Re = Ul c /i^, (4) 

3. Prandtl #, Pr = i/ k /« = ,u.C p /K (5) 

where: = }i./p = momentum diffusivity 

(kinematic viscosity) (6) 

« = K/Cp p = thermal diffusivity (7) 

u. = dynamic viscosity 

p = density 



Cp = specific heat at constant pressure 
l c = Re character i st i c length 
and in general form: 

Ru = function of Re, Pr = f (Re,Pr) (3) 

In fully developed turbulent flow systems, analysis of heat 

transfer processes is enhanced by a simple analogy between heat 

and momentum transfer which allows prediction of h from knowledge 

of shear stress (friction) acting on the fluid adjacent to the 

wall. Known as the Reynol ds-Col bur n Analogy, it is described by 

Rohsenow in the following manner C2,p.l851: 

"the heat transferred to a surface divided by the 
maximum amount which could be transferred in 
bringing the fluid to the temperature of the 
surface equals the momentum transferred in 
passing over a surface divided by the maximum 
amount of momentum which can be transferred in 
stopping the fluid relative to the surface. " 

This analogy, exactly true in parallel flow when Pr = 1 and 

reasonably accurate at Pr values not very high or very low, 
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may be used in either constant heat -flux or wall temperature 



conditions and may be expressed mathematically 


as : 


j Colburn factor = StPr^ 7 "' = C^/2 


(9) 


where: Stanton #, St = h/.pCpU 


(10) 


C f = T w /(pU 2 /2) = friction coefficient 


(11) 



T w = apparent shear stress acting on the -fluid 
at the wal 1 

Significant flow experimentation has determined that for low 
turbulent flow conditions over a smooth plate at a distance x 



from the leading edge: 




C, = .0592 Re ~ ‘ 2 

*flat plate x 


(12) 


Similarly, for turbulent flow conditions in 
pipe of diameter D: 


a cylindrical 


C, = .046 Re n - * 2 

fpipe flow u 


(13) 


Substitution of these relations into equation 
the determination of h: 

hj , flat plate of length 1: 


(9) all ows for 


StRr 2/3 = C f /2 

h w F'r 2/3 = .0592 (Re“* 2 ) /2 

p u 

h x = .0296 pC p U (Re x ) “* 2 Pr“ 2/3 

= .0296 pC U(Re w )“‘ 2 F’r -273 x K 
P K x 
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substituting -for density, the equation (6) identity, p = 



h w = .0296 K U x «£_ (Re„ ) * 2 Pr 2/3 

A r -J_j a 

K -Vy K 

recognising equations (4) and (5) identities Re v = 



Pr = /tCp /K : 



h„ = .0296 K Re.* 8 Pr 1/3 



integrating over length 1 to determine mean h^ 



Ux and 
v k 



= 1 F 1 .0296 K (U w /i/,.)* 8 Pr 1 /3 dx 

| “* y\ r*. 

1 J 0 X 



.8 1 



= .0296 K (U/v k )'° Pr 



n — 



o 



= . 0296 K (U/v,.)*° Pr 

0.8 1 



= .037 K(Re 1 )* 8 Pr 1/3 



.8 Pr 1/3 cx-B-jl 



U 



(14) 



A similar substitution of equation (13) -for turbulent -flow in 
a cylindrical pipe yields: 



h^, pipe -flow of diameter D: 



h D = .023 K(Re D )* 8 Pr 1/3 
D 



(15) 



These particular relations form the basis of many heat 
transfer correlations commonly used, rigorously applicable in 
fully developed, steady turbulent flow within established Re 
and Pr limits where temperature difference, AT, 
boundary layer is not extreme. 
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across the 



2. MODELING CONVECTIVE HEAT TRANSFER IN RECIPROCATING ENGINES 



In reci procat i ng engines, the significant process by 
which termal energy is lost is by forced convective heat 
transfer from the gas medium across cylinder inner surface 
boundaries. Large variations in the magnitude and direction 
of convective heat transfer exist from engine to engine and 
point to point within an individual engine. The conditions 
influencing heat transfer are so complex and locally 
different through an engine cycle that rigorous mathematical 
i nterpretat i on is invalid. Radical variations in 
instantaneous flow patterns allow gas velocities to be only 
vaguely character i zed . Non-uni f ormi ty of physical properties 
in the gas medium and surface materials exists. The 
combustion process is particularly difficult to model and 
further complicated by compression ignition. Two of the most 
serious obstacles question the fundamental assumption of 
fully developed, steady turbulent flow and the concept of a 
finite, instantaneous value of h. 

In operating engines, the entire cycle is completed 
within milliseconds, making accurate measurement of local 
instantaneous heat flux, q, in itself quite difficult. 
Turbulent intensity in four stroke, Cl engine cylinders is 
produced by the suction velocity of the induction stroke, the 
intensity being roughly proportional to it. To optimize 
fuel— air mi cr o-mi x i ng , off-set induction ports may be 
incorporated to create swirling gas motion. In the 
compression stroke within a simple geometry cylinder, 



induction turbulence is suppressed. 



For swirl enhanced 



cylinders with combustion chambers cut into the piston crown, 
conservation o-f angular momentum in the cup created by 
viscous drag accelerates swirl to significant angular 
velocities as the piston reaches TDC. Radial "squish" 
velocities are created near TDC and combine with swirl to 
produce a complicated toroidal motion. In the expansion 
stroke, rapid deaccel erat i on of rotational velocities in 
direct injection cylinders are superimposed by ill-defined, 
small combustion induced velocities as hot combustion product 
gases expand; in indirect injection cylinders, prechambers 
deliberately channel and exaggerate the charge velocities to 
blow out of energy cells at levels up to 500m/sec C31. In 
the exhaust stroke, the combustion products are pushed out by 
the piston. The accumulation of these factors indicate that 
turbulent boundary layers are not steady, but are 

periodically created and destroyed in a manner beyond our 
ability to model. In addition, because of the thermal 
capacity of the gases, there is a thermodynami c phase lag 
markedly apparent in reci procat i ng engines; AT and q do not 
pass through zero at the same instant indicating values of 
h = ±co and highly fluctuating values immediately about it. 

Recognizing these obstacles, i nvest i gators have 
concluded that within obtainable accuracy it is still best to 
model heat transfer as a quasi— steady process and to 
correlate cyclic observations in terms of conventional Nu, 
Re, Pr numbers. Values of the various constants and 



exponents are determined in individual engine tests; a wealth 



of experimental data provides guidance for constant (s) 

selection for operating engines or design studies of similar 

construct i on . Hohenberg states C4, pg. 27833: 

"It would be erroneous to look for a formulation 
which includes all possible influences on the 
heat transfer. This would only render the 

calculation more complicated rather than more 
accurate. The aim... must be to find out and 
mathematically describe the really essential 
factors influencing heat transfer in the engine." 

If the time for change in bulk stream conditions to diffuse 

through a boundary layer is small relative to other 

significant times in the engine cycle, then boundary layers 

may be considered quasi -steady and at any instant be that 

associated with existing conditions. Le Feuvre et al C53 

approximating diffusing time as and modeling gas motion 

in solid body rotation, successfully performed such 

calculations in a direct injection, swirl enhanced, Cl engine 

operating at high RFT1 (2000 RF'M) . As for the existence of 

phase lag, Annand C63 states that the actual quantity of heat 

transferred during this period is minute; the error arising 

from the use of a finite h at each instant of the piston 

stroke is negligible. 

The majority of the heat transferred within a 
reci procat i ng Cl engine is due to turbulent forced 
convection. On the basis of experimental data and theoretical 
consi der at i ons , this - relationship is expressed by the 

f ormul a: 



= const ant (R Q ) n F'r* 



Nu 



(16) 



Recognizing that F'r ss 0.7 -for common gas mediums due to 
predominance of 0 7 and N^, substitution of the ideal gas law 
identity, p = F'/R’Tg, into an expanded equation (16) yields: 

h = constant K>T n P n U n Tg -n l c n_1 (17) 

Assuming that actual gas velocities will be at least 
approximately proportional to mean piston speed, a -function 
o-f Sp is used with the multiplicity constant determined by 
engine design. Tg is used vice the average, (Tg + Tw)/2, as 
the calculation o-f Tg -from measured instantaneous pressure 
and gas law application already represents an average -for the 
cylinder contents. A consi der at i on o-f gas di sassoci at i on is 
neglected, as any attempt to quanti-fy it would be eclipsed by 
the previous assumption of uniform properties. During the 
combustion process, gas product composition is assumed to 
equal the burned composition. Tables and Figures of thermal 
properties (>l, Cp , F'r, 3T) for individual gases and C n H 2 n lean 
hydrocarbon fuel combustion products are provided in Appendix 
C 116,711. Cyclic local instantaneous surface temperature 
profiles, from which actual measurement of q is generally 

Y 

derived, are demonstrated by Oguri et al C83 to show 
transient magnitudes typically less than 50° F for sooty 
surfaces and less than 20°F for clean surfaces. With the 

exception of exhaust valve/port entries which may typically 
see temperatures greater than 1000°F, the temperature range 
exhibited by cylinder sleeve, head, and piston crowns are 
fairly consistent allowing assumption of a mean Tw. 



WOSCHNI 'S CORRELATION FOR MEAN INSTANTANEOUS HEAT 



TRANSFER COEFFICIENT, h„,- IN Cl ENGINES 

A recognized industry standard -for -four stroke, direct 
injection, quiescent, Cl engine heat transfer calculation is 
Woschni 's correlation, stated in U.S. units as follows C9H: 

h e = 15.480 F“ 8 CC 1 Sp+C 2 (VTg c /P c V c ) (P-Po) :* 8 B _ * 2 Tg”- 53 (18) 
where h Q = BTU/hr-f t'*— °R 

P , V c , Tg c = instantaneous values at ignition start 
Po = motoring pressure 

Based on suppressed swirl engine designs, the values of 
constants C^ and C^> are: 



i nduction 


stroke 


c i = 


CO 

i—l 

■ 

o 


c 9 = o 


( 18a) 


compression stroke 

/ 


Cl = 


n no 
sL m 


C, = 0 


( 18b) 


expansion 


stroke (motoring) 


o n 
M ►- 


— 2.2 !8 

= 0 


(18c) 



4. combustion-expansion stroke C^ = 2.28 



( 1 8d ) 

C^ = 5.9054 x 10 3 f t/sec — 0 R 
5. exhaust stroke C^ = 6.13 C 0 = 0 (18e) 

From equation (17) , a n exponent of 0.8 was determined by 
Woschni by successively changing one parameter at a time over 
maximum possible range. A consensus opinion of all recent 
published correlations is that 0.7 £ n ±= 0.8, reflecting 

agreement with the laws of similarity. Woschni developed the 
motor i ng/f i r i ng correlation and scavenging correlation on 
suppressed swirl design cylinders stating C9 , p . 3071 3 : 



"The gas motion generated during the induction 
phase dies away with time, owing to internal 
■friction. In the compression stroke, however, 
dependent upon the piston or cylinder cover 
shape, an additional gas motion can be 

induced Csqui sh— swi r 1 3 whose intensity is also 
proportional to piston speed. ' In this 

particular case, the engine involved had 

direct injection and straight -forward 

combustion chamber walls, so that decay due to 
internal -friction prevailed." 

A Tg - * 00, exponent is due to Woschni 's assumption o-f air and 

— — A — — 75 

gas product >i(Tg) « Tg and K<Tg) « Tg at n exponent = 0.8 

(see equation (17)). To approximate ill-de-fined velocity 

increase due to combustion, the difference P-Po which, unlike 

Tg-To, sees a rapid decay upon combustion completion, is 

used, multiplied by instantaneous cylinder volume, V, and 

related to the charge mass at beginning o-f combustion, 

m_ = PV/RT g _ . 
c c c 

It is at this point that Woschni makes a critical 
simplifying assumption. Recognizing the significant effect 
of luminous flame radiation heat transfer in Cl engines, 
Woschni arbitrarily accounts for combustion radiant heat 
transfer by affixing the value C 2 to equate actual measured Q 
integrated over the stroke with Q combustion determined by 
the correlation. Woschni states C9,p. 30753: 



"There is no doubt that such an influence occurs 
in the engine, and it is taken into 
consi der at i on in lumped form in the combustion 
term as defined here, but it seems impossible to 
the author to separate it from the convective 
term. The only question that remains to be 
considered, therefore, is the error made regard 
to the tendency. .. for the radiation. 
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In essence, -for application to the New Cycle Engine 
expansion cylinder operating at significantly lower peak 
pressure conditions, laws of similarity allow significant 
extrapolation of Woschni's correlation for forced convection 
heat transfer only. The correlation must be corrected to 
account for the fact that radiant heat transfer is not a 
function of <P-Po)*®. 

4. CALCULATION QF ORDER OF MAGNITUDE REDUCTION IN 

COMBUSTION FORCED CONVECTION HEAT FLUX 

The representati ve firing cycle P-Tg versus crank angle 
plot published by Woschni for natural aspirating operation is 
shown as Figure 2-1 C93. The range of additional principle 
operating parameters over which the motor i ng /f i ri ng 

correlations were developed well representive the essential 
parameters of New Cycle Engine (NCE) expansion cylinder 
operation C91: 

Parameter Woshni <W) engine NCE Expansion Cylinder 
RPM 900 max 818.3 

Sp 9 m/s max C29.6 f/s3 33.33 ft/sec 

Supercharge 3 max 1.0 

pressure 
rat i o 

Tw nominal 500° K C900°R3 900° R 

L,B unknown 1.22204 ft 

Bore and stroke dimensions were not provided, however from 
known maximum data point values Sp = 9 m/sec and = 900, 



the approximate stroke length L may be determined: 



h g proportional to B ’ ^ in forced convection heat 



L x 60 Sp x (60) (9) x 0.3 m CO. 98 ft: 

2Npp|y| (2) (900) approx. 

For conventional Cl engine geometry, bore diameter B is 

approximately equal to L; therefore for effective comparison 

of 

transfer, consi derat i on of bore diameters of similar 
magnitude dimensions may be neglected. The significant 
difference between combustion-expansion stroke profiles are; 

1. An order of 15 magnitude increase in peak 
combustion pressure in Woschni's 
conventional high r , Cl engine. 

2. A reduced Tg decay rate, thus increased 
(Tg-Tw) forced convective heat transfer 
potential in the New Cycle Engine 
expansion cylinder. 

3. A delayed fuel i n j ec t i on/ i gni t i on start 
in the New Cycle Engine expansion 
cyl i nder . 

Assuming that the magnitude of mean instantaneous gas 
velocities through the expansion strokes are roughly equal, 
the approximate ratio of forced convective heat transfer loss 
in the respective engines may by obtained by determining the 



mean ‘ratio of instantaneous heat flux, 



over a 90° crank 



angle domain from ignition start. A selection of 90° is 
based on; 



1. Establishing a normalized basis in which to 



compare the order of magnitude reduction in 



forced convection versus radiant heat 

transfer; the crank angle domain over which 
luminous flame radiation is significant in 
medium speed Cl engines less than 1000 RF'M 
is nominally 90° from combustion start CS]. 

2. At 90° the slope of pressure and 

temperature decay transients induced by 
combustion in both engines become small. 

Setting the radiant heat transfer constant C, ? = 0 in equation 
( 18 d ) : 



q W /q NCE 



q W /q NCE 



« h e <Tg-Tw) w /h e (Tg-Tw> NCE 

Up. 8§p- S /B* 2 Tg* 53 II CTg-Tw3 w 

« — 

CP* 8 Sp* 8 /B* 2 Tg * 53 : CTg-Tw3 NCE 



Neglecting B consideration and assuming equal magnitude gas 
velocities: 

q w c p. 8 /Tg.53 ][Tg _T w ] w 

« (19) 

^NCE CP- S /Tg- 53 i;Tg-Twn NCE 

where Tw = 900° R 

Table 2-1 summarizes the calculations based on Figure 2-1 and 
Appendix I data at 5° increments from combustion start. The 
mean heat flux reduction ratio may then be determined by 
numerical integration using Simpson's rule: 
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q w - _1_ Cq w /q NCE JdU (^.U) 

q wr r A90° J 

'^■forced convection u 



q w 


^ 5 ° q w q w q w q w q w q w 

Cl +4 +2 +...+2 +4 +1 3 

3 q NCE q NCE q NCE q NCE q NCE q NCE 

0° 5° 10° 80 85 90° 


q NCE 


A90° 


= 


A5°/3 C 328 . 67 3 = 6.087 approx . 

90° 






TABLE 2-1 



Calculation of order of magnitude reduction in combustion 
forced convection q from Woshni's test engine to New Cycle 
Engine expansion cylinder operation at 5° crank angle 
i ncrements 



Crank angle 




Noshni 


engine (N) 




NCE Expansion Cylinder 


V q NCE 


past ignition 


p 2 
kp/cr 


Tg 


P 


Tg 


P 


Tg 


equation (19) 


start 


•K 


ate 


•R 


ate 


♦R 




0 


40 


875 


38.71 


1575.0 


4.15 


1896.6 


4.46 


5 


60 


1400 


58.07 


2520.0 


4.25 


2252.6 


9.14 


10 


63 


1525 


60.97 


2745.0 


4.09 


2571.2 


9.29 


IS 


60 


1600 


58.07 


2880.0 


3.92 


2864. 1 


8.68 


20 


55 


1650 


53.24 


2970.0 


3.76 


3136.0 


7.94 


25 


47 


1630 


45.49 


2934.0 


3.26 


3064.2 


7.92 


30 


40 


1600 


38.71 


2880.0 


2.89 


2998.7 


7.68 


35 


33 


1540 


31.94 


2772.0 


2.59 


2939.1 


6.66 


40 


26 


1480 


25.16 


2664.0 


2.34 


2884.8 


6.2 


45 


22.5 


1440 


21.18 


2592.0 


2.13 


2835.2 


5.76 


50 


20 


1400 


19.36 


2520.0 


1.96 


2790.0 


5.65 


55 


17 


1365 


16.45 


2430.0 


1.81 


2748.8 


5.16 


60 


IS 


1330 


14.52 


2394.0 


1.69 


2711.3 


4.92 


65 


12.5 


1305 


12.1 


2349.0 


1.58 


2677.1 


4.46 


70 


11 


1280 


10.65 


2304.0 


1.49 


2646.0 


4.17 


75 


10 


1255 


9.68 


2259.0 


1.42 


2617.8 


3.97 


80 


9 


1230 


8.71 


2214.0 


1.35 


2592.3 


3.75 


85 


8.5 


1215 


8.23 


2187.0 


1.3 


2569.2 


3.68 


90 


8 


1200 


7.74 


2160.0 


1.25 


2548.4 


3.59 




FIGURE 2-1: Woschni's test engine firing cycle; natural 
aspiration [9] 
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ESTIMATION OF MEAN GAS VELOCITY MAGNITUDE, 0 
5. 1 CONSTRUCTION OF CHARGE IMPULSE VELOCITY , Ui , TRANSIENT 

Acceptance of the general applicability of Woschni 's 
correlation to the New Cycle Engine expansion cylinder 
balances on the assumption of comparable gas velocity 
magnitudes. While geometric flow patterns are indeterminate 
in preliminary design, the magnitude of the charge impulse 
into the expansion cylinder leaving the regenerator port can 
be ascertained. Recognising the inevitable deacceleration of 
this impulse in a simple geometry cylinder on the expansion 
stroke, a conservative estimation of the mean bulk gas 
velocity at the walls, Uexp, through the expansion stroke can 
be made with minimum error. Substituting Oexp directly into 
the steady, fully developed turbulent flow equations (14) and 
(15) for flat plate and pipe flow should, if the magnitudes 

of cylinder velocities are in fact similar, be reasonably 

• 

close to h e calculated from Woschni 's correlation. To reduce 
uncertainties, the rel at i onshi ps will be first compared in 
motoring conditions from TDC to 45°, prior to ignition start. 

Dent et al C 1 0 1 successfully used a similar method in 
approx i mat i ng forced convective heat transfer loss in a four 
stroke, direct injection, high swirl, Cl engine operating at 
peak TDC swirl ratios equal to 20. Predictions using 
Woschni 's motoring correlations failed uniformly low compared 
to measured, spatially averaged, instantaneous heat flux. 
This is expected, considering the central assumption of 
suppressed swirl inherent in application of Woschni 's 
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correl at i on . 



By substituting directly into the flat plate h 



equation (14), measured instantaneous angular gas velocity u> r 
at l c equal to bore circumference (2irr), remarkable agreement 
was obtained. Agreement was not as good at firing conditions 
due to the transient combustion processes. 

From ideal gas law, at a given angle 6 in the expansion 
stroke, charge mass, m, in the expansion cylinder: 

PV (2116.224 lbf /f t 2 -atm> PV 

m 9 lbm = ~~T = = 

RTg ( 53 . 90 ) Tg 

= 39.262 PV (21) 

Tg 

Taking 5° crank angle as the unit increment of time at 818.3 
RPM; the charge mass rate, entering the expansion cylinder 
is: 

Am (5° ) 

♦ - - 

m<j i l _ as Am (5°) =15 
A lbm/sec 

At (5°) (N RpM /60) 360° 



Am (5° ) 

.0010184 sec 



( 22 ) 



The continuity equation for impulse flow leaving the 
regenerator port cross-section area, Ap , equal to 0.80 ft^ - 
is: 



m 0 



= p U i A 



P 



(23) 



Therefore, by rearranging equation (23), impulse velocity 
leaving the regenerator port, Hi, can be determined: 



m 0 (53. 90) Tg 

(21 16.224 1 b-f /+' t 2 -atm) P ( . 080) 



.31837 m 0 Tg 



m 0 RTg 



Ui 



•f t/sec 



PA. 



A summary o-f the calculations using instantaneous P, Tg, V 
values -from Appendix A at 5° increments is provided in Table 
2-2. Charge mass as a -function o-f expansion stroke crank 
angle is plotted in Figure 2-2, indicating that the main 
impulse -from the compression cylinder via the regenerator 
exit port completes 10° past ignition start. 

At 45° combustion start 7.m. . = 2. 1905x 10 -i = 89.37. 

total _o 



2.4537x10 ^ 



At 



55° 



main impulse completion 

7. m 



total 



2. 28394 x 10_ J 
2. 4537 x 10 " 2 



rn “ r TV 

— *7 • J. . O /» 



Ui as a -function o-f expansion stroke crank 
in Figure 2-3, indicating a peak Ui equal 
near -fuel injection start (30°). 



angle is plotted 
to 519.95 -ft/sec 
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TABLE 2-2 



Calculation of charge Hi leaving the regenerator port into 
the New Cycle Engine expansion cylinder at 5° crank angle 
i ncrements 



Crank 






Angl e 


m 




m 


Ui 


e° 


lbm 


AlbmE-3 


1 bm/sec 


ft /sec 




eq. (21) 




eq. (22) 


eq. (24) 


TDC 


0 


0 


0 


0 


5 


1 . 5550E-4 


0 . 1 555 


0. 1527 


41.82 


10 


7. 7090E-4 


0.6154 


0 . 6043 


133. 14 


15 


2 . 1 309E-3 


1 . 360 


1 . 3275 


236 . 92 


20 


4. 5821E-3 


2.451 


2 . 4069 


352.69 


25 


8. 31 13E-3 


•—> a 7 L 


3. 6618 


458. 13 


30 


1 . 2936E-2 


4.625 


4.5413 


519.95 


35 


1 . 7373E-2 


4.437 


4.357 


499. 15 


40 


2. 049 IE— 2 


3.118 


3. 0614 


382.68 


45 


2. 1905E-2 


1 .414 


1.3882 


202.03 


50 


2. 2887E-2 


0 .9816 


0.9641 


162.73 


55 


2. 2839E-2 


0. 0028 


0 . 00275 


0 . 55 


60 


2. 2932E-2- 


0.0426 


0 . 04 1 83 


9. 72 


65 


2. 2992E-2 


0. 0596 


0.05852 


15.54 


70 


2.3056E-2 


0. 0641 


0.06294 


18.83 


75 


2. 3307E-2 


0 . 25 1 0 


0.2465 


81.44 


80 


2. 3512E-2 


0.2049 


0 . 20 1 2 


72. 76 


85 


2. 3680E-2 


0. 1687 


0. 1657 


65. 12 


90 


2. 3820E-2 


0. 1397 


0. 1372 


58.15 


95 


2. 3937E-2 


0. 1167 


0.11 46 


52. 03 


100 


2 . 4035E-2 


0.0981 


0 . 09633 


46.56 


105 


2. 41 18E-2 


0.0827 


0.08121 


41.54 


110 


2.4188E-2 


0 . 0702 


0. 06893 


37. 12 


115 


2. 4247E-2 


0. 0596 


0.05852 


33 . 0 


120 


2. 4298E-2 


0 . 05 1 1 


0 . 050 1 8 


29.84 


125 


2. 4342E-2 


0. 0436 


0.04281 


26. 13 


130 


2. 4379E-2 


0.0372 


0. 03653 


23.04 


135 


2.441 IE-2 


0 . 0320 


0 . 03 1 42 


20.41 


1 40 


2.4438E-2 


0 . 0270 


0 . 0265 1 


17.67 


145 


2. 4462E-2 


0.0239 


0. 02347 


15.99 


150 


2. 4480E-2 


0 .0183 


0. 01797 


12.48 


155 


2. 4496E-2 


0 . 0 1 58 


0.01552 


10.94 


160 


2. 4509E-2 


0 . 0 1 30 


0. 01276 


9. 12 


165 


2.4519E-2 


0 . cl) 1 0 1 


0. 00992 


7. 15 


170 


2. 4527E-2 


0. 0074 


0 . 00727 


tr nn 

U m AD 


175 


2. 4532E-2 


(D . 005 1 


0 . 0050 1 


3.65 


BDC 


2.4537E-2 


0.0052 


0 . 005 1 1 


3. 70 



45 



I bra 



.025 



. 020 -- 



. 015 - 



. 010 - 



. 005 -- 



0 




TDC 



0 ! °c|e9rees 



FIGURE 2-2: 



Expansion cylinder mass charge, m, 
expansion stroke crank angle 



versus 
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ft /sec 



600 



400 - 






200 — 



0 




TDC 



0 ° decrees. 



FIGURE 2-3: Expansion cylinder charge\/i versus expansion 
stroke crank angle 
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VERIFICATION OF TURBULENT FLOW 



To verify 


turbul ent 


flow conditions 


I eavi ng 


the 


regenerator port , 


cal cul ated 


Re at 30° using l c 


equal to 


the 


regenerator port 


cross-sect i on minor dimension 


(depth ) 


of 


0.75 inch is: 










Re 30° 


Ui 1 c Ui 1 c 

Vi, JJ./p 


• Ui l c P 

^ RTg 


(25) 





From Appendix C, Table C-l , ^(Tg) and Cp(T) for air is: 

/i lbm/ft-sec = 2. 163x 10 -7 (Tg ) ' 645 (26) 

(Cl ) 

C BTU/1 bm-hr = 0.224 + 0. 262x1 0” 4 Tg (27) 

P ( C2 ) 

Substituting instantaneous values of P, Tg -from Appendix A 
and >t(Tg) identity equation (26) into equation (25): 



P -^-,0 = 6.10171 atm 
Tg 30 o = 2194. 29°R 




(519.95) (0.75/12) (2116.224 1 bf /f t^-atm ) (6.10171) 
(2. 163xl0 -7 ) (2194.29) ' 645 (53.90) (2194.29) 



= 1 



1 476x 10 5 



turbul en t 
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5.3 ESTIMATION QF EXPANSION STROKE, Uexp 

The calculation of the mean impulse velocity, Ui , -for 
the purpose o-f Uexp approxi mati on , is correctly restricted to 
the main impulse domain TDC to 55° only, neglecting residual 
107. charge mass regenerator bleed. Combustion induced 
perturbations during crank angles 45°-55° may be neglected. 
Ui may then be calculated using numerical integration by 
trapezoidal rule and is equal to: 

Oi r f / _ _ _ = 1 f 55 ° Ui d0 (28) 

TDC— 55° A55°J tdc 

(A5°/2)C Ui TDC +2Ui 50 +2Ui 1C)0 . . .+2Ui=, 00 +Ui 55 o 3 

A55° 

% 

= ( A5° /2 ) C 5979 . 03 3 = 271.77 -ft/sec 

A55° 

A reasonable assumption must be made as to the expected 
value o-f Uexp over the entire stroke. The principle -factors 
causing charge velocity deaccel er at i on include cross-sect i on 
area expansion exiting the regenerator port, internal 
-friction, cylinder volumetric expansion, and angular momentum 
deaccel er at i on o-f any existing rotating -flow patterns. 
References C5, 10, 11, 12, 133 unanimously agree that in swirl 

enhanced engines, peak instantaneous swirl at TDC is much 
greater than the mean swirl valves over the stroke. 

Le Feuvre et al C53 calculated -for a 4.5 inch bore direct 
injection, Cl engine operating at 2000 RF'M, that 107. velocity 
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decrease occurred each stroke due to viscous dissipation 
alone. A conservative estimate of Uexp = (Oi 700 - 550 > /2 is 

therefore assumed; selection of this value should reduce any 
error in either direction. It is noted that ^joc-BDC ^ rom 
Table 2-2 in firing conditions equals 102.9 ft/sec. 

Combustion induced velocities due to charge gas expansion in 
spark ignition (SI) engines have been calculated by Annand 
C3J to be on the order of 33 ft /sec, with expected velocities 
in Cl engines less due to diffusion flame conf i gurat i on . 
Thus, let Uexp = (Ui TDC _ 550 ) /2 = 135.9 ft /sec. 

6 . COMPARISON OF h c CALCULATED FROM NQSHNI'S CORRELATION 

AND THEORETICAL MODELS 
6 . 1 EXPANSION STROKE; MOTORING CONDITION 

Using the above first appr ox i mat i on of Oexp, estimates 
of h e from theoretical equations (14) and (15) for steady, 
fully developed flow may be reached. Expanding pipe flow 
equation (15) at 1 = B and U = Uexp yields: 

h D = . 023 (K/B) Re D * 8 Pr 1/3 = . 023 (K/B) CUexpB/v k 3 * 8 Pr 1 /3 

Substituting the identities = /t/p , K = r, and 

p = p/RTg: 

h D = . 023C py H.CPUexp/ y wRTg3 ■ S B _ - 2 Pr” 2/3 

Substituting the C (T) and >i(Tg) identities from equations 
(27) and (28) and appropriate conversions: 



h D ” 



. 023 C . 224+ . 262 x 10 4 Tg3C2. 163x 10 7 (Tg) * 645 <3600sec/hr ) 3* 
[<2116.224 lbf/ft 2 -atm P Uexp) / <2. 163x 10 _7 Tg‘ 645 R?g) 3 * 8 * 
B ,A Pr ^ / - 



which simplifies at known R, Dexp, Pr values to: 



> \ ■ 8 



*D 



BTU/hr-f t i -°R 



= 9. 194C2. 24+2. 62x 10 4 Tg3 



<135.9 P) 



B* 2 Tg* 671 



(29) 



An identical expansion of the flat plate equation <14) at 
1 = irB , expansion cylinder ci rcumf erence , yields: 



= 11.674C2. 24+2. 62x10 4 Tg3 



<135. 9 P) * 13 



BTU/hr-f t^- 0 R 



B*^ Tq 



. 671 



(30) 



A summary of h calculations for the charge induction 
monitoring angles, TDC-45 0 , comparing equations <29) and (30) 
to Woshni's motoring correlation equation (18c), is provided 



in Table 2-‘ 



The agreement is quite good, even with 



acceptance of reasonable error in Uexp approximation. 
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TABLE 



Comparison of R e calculated by Woshni's correlation and 
theoretical models during charge induction into the New Cycle 
Engi ne. expansi on cylinder 



Crank 



Angle 


h e 


h D 


h l 


e° 


eq . ( 18c ) 


eq. (29) 


eq. <30 


5 


16.994 


15.221 


19.477 


10 


20. 191 


18. 143 


23. 215 


15 


24. 006 


21 . 593 


27 . 63 1 


20 


28.209 


25 . 396 


32.496 


25 


32. 029 


28. 839 


36.902 


30 


34.236 


O • 79 o* 


39. 403 


35 


33. S88 


30 .414 


38.918 


40 


31 . 168 


27.897 


35 . 697 


45 


27.164 


24.251 


31 . 037 



h in BTU/hr— f t 0 R 



6.2 



EXHAUST STROKE 



Attempts to directly estimate the mean bulk gas 
velocity at the walls during the exhaust stroke, Uexh, in 
light o-f preliminary design restrictions and post-combustion 
ambiguities, are not valid. 

An indirect method allowing h calculation is to correct 
equation (29) by a -factor equivalent to the 

scavenging/motoring Uexh correction observed by Woschni -for 
suppressed swirl geometry, direct injection, Cl engine 
cylinders. From equations 18a and 18c: 

Uexh correction -factor = (6.18 Sp/2.28 Sp ) * 8 - 2.22 

When the regenerator exit valve 3 opens, the 

signi-ficant -flow pattern will be parallel to the bore axis, 
up the sleeve, and through the regenerator port; -flow 
direction is therefore analogous to pipe -flow. From Appendix 
C, the identities /(.(Tg), Cp(Tg), Pr (<f) , and H' < T g ) -for lean 
C | _ l H- :>n hydrocarbon -fuel combustion products at P = 1 atm, 

$ = 0.75 may be calculated: 

1 . 47 x 10 — 7 (Tg ) m7 1 . 47 x 1 0 — 7 ( Tg ) ml 

/(. lbm/ft-sec = = 

1 + .027$ 1 + .027 (.75) 



= 1.441x10 7 (Tg ) ’ 7 (31) 

(C3) 

C BTU/1 bm— °R = 0.247 + 0. 243x1 0 -4 Tg (32) 

P (C5) 

Pr(3T) = 0.05 + 4.2<y-l) - 6. 7(^-1) 2 (33) 

( C4 ) 

Hy'(Tg) s Appendix C, Figure C-l (34) 



Applying the Uexh correction factor = 2.22 and the identities 
/UTg) , Cp(Tg) defined above to pipe flow equation (29) 
yi elds: 

h D = 2. 22 (.023) ( . 247+.243X 10 -4 Tg) < 1 . 441 x 10~ 7 Tg‘ 7 ) <3600sec/hr ) * 
C <2116.224 lbf/ft 2 -atm P Uex p ) / ( 1 . 441 x 10 _7 Tg * 7 RTg ) 3 ‘ 8 * 

B - Pr 

which simplifies at known R C p, Uexp values to: 

135.9 P 

h n = 14.835C2.47 + 2.43x10 4 Tg3 Pr 2/3 

BTU/hr-f t A -°R B'^Tg* toto 

(35) 

From Appendix A; P = 1.05 atm and 2344. 58°R ^ Tg i 2439. 06° R 

TDC BDC 

through the exhaust stroke. A comparison of h calculated 
from equation (35) and Woschni's exhaust stroke correlation 
equation (18e) follows: 



crank angle 


Tg 

°R (°K) 


«T ( T g ) 
eq. (34) 


Pr (j*> 
eq . ( 33 ) 


h e 

eq. ( 18e) 


h D 

eq. (35) 




BDC 


2439.06 

(1355) 


1.275 


0 . 698 


17.587 


17.062 




TDC 


2344. 58 
( 1303) 


1 . 280 


0. 701 


17. 959 


17.331 










h 


in BTU/hr- 


f t 2 -°R 




The above 


series of 


calculations suggest that 


to a first 


approximation, application 


of Woschni's correlation 


i n 



determining forced convective heat transfer loss for the New 
Cycle Engine expansion is valid. 
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CHAPTER THREE 



CALCULATION OF UNCQRRECTED HEAT TRANSFER 
LQSS/HEAT RELEASE PER CYCLE, ft 

Table 3.1 summarizes the calculation of total 

instantaneous heat transfer rate, Q, through a complete 360° 

firing cycle of the New Cycle Engine expansion cylinder using 

Woschni's correlation for hi . For calculation of FL in 

e g 

combustion, motoring F* is generated graphically in Appendix 
A, Figure A-l: 

QgTU/sec = *^e A HT *Tg-7w) (1 hr/3600 sec) (36) 

Due to the establishment in preliminary design of a simple 
cylinder geometry, A^j as a function of crank angle may be 
determined directly from instantaneous V values; defining 
instantaneous gas volume height as L^ and neglecting cylinder 
sleeve clearance: 



which 



V = HB 2 
4 



II 

H 

ZC 

<r 


2ttB j! ' + ttBL 0 = 
4 


2ttB^ + 
4 


at B 


= 1.22204 


ft: 




a ht 


= 2.346 


+ 3. 


273 C V 3 



irDC4V/irB 2 3 



(37) 



The distribution of h 0 and Q as a function of crank angle is 
plotted in Figure 3—1. It is immediately clear studying the 
h e distribution over the combustion stroke 55° to BDC , that 
Woschni's correlation fails to correctly extrapolate luminous 



5S 



f 1 ame 



operat i ng 



radiant heat transfer loss to low r c 
conditions. No increase in h 0 occurs due to the predominant 
influence of increased temperature: associated gas 1 ayer p 

decrease and increase overcome the effect of K increase, 
thus lowering h g , modeled by laws of similarity solely on 
forced convective heat transfer. The low r c and delayed fuel 
injection operating conditions effectively eliminate any 
remarkable pressure rise in the cylinder, already at 45° 
down-stroke at ignition start. Recalling Woschni's coupling 

Q 

of radiant heat transfer in combustion to C/ ? <P-P ) , the 

desired "lumped form" term effectively disappears. 

Uncorrected total heat transfer loss per cycle, Qy, is 
determined by numerical integration under the 0. curve using 
Simpson's rule: 



Q 



T 



BTU/cycl e 



* * 

} BDC + 4Q 185' 



= . 00 101 84 [0^+ 4Q , aero + 2Q 



1 90* 



+ 20 . -,,-, 0+40 



170° 175° BDC 



+ 0 , 



(38) 



= . 0Q 1 0 1 84 C 7 1 4 1 .6133 = 2.42434 BTU/cycle 



From Table 2-2, total charge mass per cycle equals .024537 
lbm. At known specific heat release per cycle, A heat, equal 
to 500 BTU/lbm charge, the ratio of total heat transfer 
loss/heat release per cycle, Tty, equals: 

Tty = 2.42434 BTU 

uncorrected (500BTU/1 bm) (. 024537 lbm/cycle) 

= 19.761*/. 
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The inability of Woshni's correlation to account for radiant 

heat transfer in low r c operation is clearly demonstrated by 

repeating the numerical integration with the combustion term 

in equation (18d) eliminated (C 7 = 0) s 

Simpson's Multiplier = 7038.088 

Q t = 2.3892 BTU/cycl e 

1l T = 19.4747. 

uncorrected 

Considering the accuracy of the correlation in general, for 
preliminary design prediction purposes setting = 0 is 

considered valid. 

To correct for luminous flame radiation heat transfer 

loss, a value of Cj ^ ^ ^ radi at i on must be determined which, 

in comparison of q^/q^^- ^ Qrcecl convection calculated in 

Chapter Two, will allow first approximation of the required 

Tl-|- correction factor. As there exists no theoretical basis 

for coupling magnitude of luminous flame radiation to 
8 

(P-P ) ' , an attempt to determine the desired ratio from 

8 

(P— P ) ' profiles must be disregarded in favor of calculation 
founded on known principles of radiation heat transfer in 



diffusion flame combustion. 



TABLE 



1 



Calculation o-f New Cycle Engine expansion cylinder -firing 

» 

cycle Q using Woschni's correlation at 5° crank angle 
i ncrements 

crank angle h e A|_|y Tg-Tw Q 

0° BTU -ft 2 0 R BTU 

o sec 

hr— ft 2 -°R 

eq(18) eq<37) Tw=900°R eq (36) 

charge induction equation <18c) 


5 


16.944 


2.357 


1249.04 


13. 856 


10 


20. 191 


2.390 


1259.25 


16.878 


16 


24.006 


•2.445 


1295.73 


21.128 


20 


28.209 


O cr OO 
■ \-J xL xL 


1325.26 


26. 190 


25 


32.029 


2. 618 


1329.71 


30. 972 


30 


34.236 


2. 734 


1294.29 


33. 652 


35 


33.888 


2.867 


1216.61 


32.834 


40 


31. 168 


3.017 


1110.57 


29 . 009 


45 


27. 164 


3.181 


996. 6 


23.924 


combust i on -ex pan si on 


equat i on < 1 8d ) 






50 


26.086 


3. 358 


1351 . 6 


32 . 888 


55 


23.969 


3.548 


1671.21 


39.445 


60 


22. 066 


3.741 


1964. 13 


45.038 


65 


20 • 85o. 


3. 944 


2236.02 


51 . 083 


70 


18.710 


4. 152 


2164. 15 


46. 700 



TABLE 3-1 (CONT'D) 



angl e 
0 


h e 


a ht 


Tg-Tw 


Q 


e 


BTU 

hr-f 1 2 -° R 


ft - * 


°R 


BTU 

sec 


75 


17.084 


4.362 


2098.73 


43. 444 


80 


15.742 


4.573 


2039. 1 1 


40.775 


85 


14. 632 


4.783 


1984.75 


38. 584 


90 


13. 715 


4. 990 


1935. 18 


36.789 


95 


12. 841 


5. 192 


1889.98 


35. 002 


100 


12. 099 


5. 388 


1848. 18 


33. 473 


105 


1 1 . 457 


5. 577 


1811.25 


32. 148 


1 10 


10. 865 


5. 757 


1777.09 


30.877 


1 15 


10.438 


5. 927 


1746. 03 


30. 006 


120 


1 0 .015 


6. 087 


1717.83 


29.089 


125 


9. 652 


6. 236 


1692.29 


28 . 36 1 


1 30 


9.329 


6. 373 


1669.22 


27.567 


135 


9. 052 


6. 498 


1648.44 


26. 934 


140 


8 • 8o>0 


6.61 1 


1629.82 


26. 428 


145 


8.632 


6. 71 1 


1613.21 


25. 959 


150 


8.467 


6.797 


1598.51 


25.554 


155 


8.334 


6.871 


1585.62 


25.221 


1 60 


8.229 


6 . 930 


1574. 46 


24.941 


165 


8. 151 


6. 978 


1564. 96 


24.725 


170 


8 . 086 


7.011 


1557.07 


24.517 


175 


8. 044 


7 . 03 1 


1550.74 


24. 363 


BDC 


8. 028 


7.037 


1545. 46 


24.260 
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TABLE 



1 (CONT'D) 



angle 


R e 


a ht 


r .: n rtp 

1 g— 1 w 


Q 


0 

0 


BTU 


■ft*- 


0 R 


BTU 




hr~f t 2 - 0 R 






sec 


haust 


equation <18e) 








BDC 


17. 587 


7. 037 


1539. 06 


52.909 


185 


17.604 


7.031 


1534.63 


52.763 


190 


17.619 


7.011 


1530. 70 


52 . 523 


195 


17.634 


6.978 


1526. 81 


52. 137 


200 


17.649 


6.930 


1522. 95 


51 . 741 


205 


17.664 


6.870 


1519. 13 


5 1 . 208 


210 


17.678 


6.797 


1515. 36 


50 . 578 


215 


17.693 


6.710 


1511.61 


49. 894 


220 


17.707 


6.611 


1507.99 


49.035 


->25 


17.721 


6.498 


1504.41 


48. 121 


230 


17.736 


6 . 373 


1500. 90 


47.125 


235 


17.748 


6.236 


1497.47 


46.037 


240 


17.761 


6. 087 


1497.47 


44. 870 


245 


17.774 


5.927 


1490.89 


43. 628 


250 


17.786 


5.756 


1487.75 


42.308 


255 


17.798 


5.576 


1484.71 


40.929 


260 


17.810 


5.388 


1481 .71 


39.498 


265 


17. 821 


5. 192 


1478. 95 


38. 012 


270 


17.832 


4. 990 


1476. 24 


36. 488 


275 


17.842 


4. 783 


1473. 65 


34. 933 


280 


17.852 


4.573 


1471 . 18 


"7* “7* T / n 

•JO a JOjL. 
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TABLE 



1 (CONT'D) 



angle 


R e 


a ht 


Tg-Tw 


Q 


0 

0 


BTU 


ft 2 


°R 


BTU 




hr — f t 2 - 0 R 






sec 


285 


17.861 


4.362 


1468.82 


31 . 786 


290 


17.870 


4.152 


1466.58 


30.226 


295 


17.879 


3. 944 


1464.46 


28.685 


300 


17.887 


3.741 


1462.44 


27. 198 


305 


17.895 


3.545 


1460.54 


25.737 


310 


17. 902 


3 . 356 


1458.73 


24.346 


315 


17.909 


3. 181 


1457.03 


23 . 057 


320 


17.915 


3. 017 


1455.41 


21 . 851 


325 


17.921 


2. 867 


1453.88 


20. 750 


330 


17.927 


2. 734 


1452.41 


19.774 


ttct 

OOJ 


17.933 


2. 618 


1451.01 


18. 923 


340 


17.938 


n icr ‘ i f ~\ 


1449.67 


18.217 


345 


17.944 


2.445 


1448.36 


17.651 


350 


17.949 


2.390 


1447.09 


17.244 


355 


17.954 


2.357 


1445. 83 


16. 996 


TDC 


17.959 


2.346 


1444.58 


16.906 
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h e : BTu/hr-ft-^ (Q : BTu/ sec) 




EXHAUST 



_ o EXPANSION 

C7 • de0rees 



FIGURE 3-1 : 



Expansion cylinder firing cycle h e and Q versus 
crank angle; Woschni’s coorelation 
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CHAPTER FOUR 



RADIANT HEAT TRANSFER IN THE NEW CYCLE ENGINE 

EXPANSION CYLINDER 

1. THEORETICAL BACKGROUND 

When a temperature difference exists between two 
bodies, part of the thermal energy is transferred as 
el ectromagneti c waves of radiation at varying wavelengths. A, 
character i st i c of the body temperature and identity. In the 
el ectromagnet i c spectrum, thermal energy is propagated 
through .01 >im *= A ^ 100 ,um. Bodies above absolute zero 
temperature continuously emit energy which, when striking 
another body, may be partly absorbed , fo, reflected, ^ , or 
transmitted , T, where: 

fr + t, + 7 = 1 (39) 

Thermal radiation reflected from and transmitted through a 
body is again absorbed, reflected by, or transmitted through 
any additional bodies ultimately establishing an equilibrium 
state. An ideal body that absorbs all radiation is known as 
a black body; if no radiation is transmitted, as in the case 
for most industrial liquids and materials, the body is known 
as opaque. Combustion chambers covered with a thin film of 
oil contaminated with soot particles will absorb a 
significant percentage of incident radiation, T = 0. 

Thermal energy emitted from a body unit surface per 



unit time at a particular wavelength, 



the monochromatic 



emissive power, E^ , can be determined as a -function o-f 
temperature T(°K) by Planck's Law: 



= 2itcc X J Ccc_ /e c b^-l 3 1 
* . 3 P P 

erg/sec-cm 



(40) 



where: 



: speed o-f light 2. 997902 x ID 10 cm/sec 

= Planck's constant 6. 623777x erg-sec 
= Boltzmann's constant 1.38026x10"*^ erg/°K 



X = wavelength in cm 
T = temperature in °K 

will increase with increasing A to a maximum value then 
decreases, the distribution shi-fting toward the shorter 
wavelengths at higher body temperatures. The relationship 
between T(°K) and X v is determined by Wein's Law o-f 
Di spl acement : 



X m _ w T(°K) = 2898 cm— 0 K 
m q a 

The total radiant emissive power, E, may therefore be 
calculated at a given body temperature by integrating the 
area under the curve: 

E = f“ E. dN (41) 

Jo X 

For a given black body, substitution o-f equation (40) into 
equation (41) yields the Stephan Boltzmann Law: 

fc Black Bod>; "o’ (°R> 

BTU/hr— ft^ 

where: cr Q = Stephan Boltzmann constant 1.713x10 ^BTU/hr— f t '~° 



64 



Real bodies, however, are not absolutely black. Those 
bodies o-f a continuous radiating electromagnetic spectrum 
where the curve is proportional to that o-f an absolute 

black body at each 7 \ are known as grey bodies. The ratio of 

emissive power of a grey body to an absolute black body is 
known as the coefficient of emissivity, £: 

E grey body (T) = £ (43) 

E black body (T) 



From an analysis of equilibrium emittance and absorptance 
between a grey body and its surroundings at the trivial case 
of each equal to the same temperature, it can be demonstrated 
that £ is equal to Jo (Kirchoff's Law). This important 
phenomenon, strictly correct only within the defined 
limitations, is fortunately observed to be approximately 
correct in many real conditions where temperature differences 
are not severe. Additionally, many real industrial materials 
radiate near grey at very high values of Jo; instantaneous 
radiant heat transfer rate, Qp , between two near-grev bodies 
1 and 2 of different temperatures at a distance apart may 
therefore be calculated from equation (44): 



°R1 



Btu/sec 



4_ t 4- 

1 ,2^o r 1 ^"HT' 1 1 lo 



= £ , oAuT(Ti* r ~Tn* T ) (1 hr/3600 sec) 



(44) 



where: £ 



1 , 



= 1 + 1-1 



1,2 



= geometric view factor between bodies 1 and 2 



65 



Tables providing formulas for £ and F calculation 
between multiple bodies at various common geometries are 
published in standard heat transfer texts. For gas or 
incandescent solid particles radiating as a defined unit body 
within an encompassing simple geometry combustion cylinder, F 
may be assumed equal to unity. 

2. NQN-LUM1NQUS GAS RADIATION 

2. 1 CHARACTERISTIC HEAT TRANSFER EQUATIONS 

Unlike solid bodies which radiate across continuous 
wavelength spectrums, gases emit and absorb energy at direct 
wavelengths only in line spectra. For simple monatomic and 
diatomic gases such as 0^,N^, and H^, the line spectra are 
negligible and may be considered transparent <T = 1). For 
complex gases such as generated in exhaust products, however, 
radiant capacity may be significant. CCU and H^O vapor 
predominate as radiant emitting gases in combustion products; 
the significant bands of radiation are C143: 

CCU H^O 

■ ■■ jL ■ xL. - 

Band 1 2.36 - 3.02 ,nm 2.24 - 3.27 /im 
Band 2 4.01 - 4.80 Rm 4.8 - 8.5 /im 
Band 3 12.5 - 16.5 jun 12 - 25 un 

In radiating gas, both emittance and absorptance exist 
throughout the gas volume. Penetrating radiant energy will be 
reduced by absorption proportional to the number of molecules 
in the ray path as a direct function of individual gas 
partial pressure, temperature, and path length, X. This 
intensity attenuation is expressed by Beers Law: 
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(45) 



t0 A = l~exp<k x X) 

where: k. = monochromatic absorption coe-ff icient per 

unit length 

X = radiant path length 

A gas volume is, of course, not grey. For ease of 

calculation an equivalent gas emissivity, £g, is defined in 
accordance with equation (46): 

£g = 1 - exp <k g l R ) (46) 

where: k g = apparent gas absorption coefficient per unit 

1 ength 

1 R = effective mean beam length 

For a defined gas volume of V and boundary surface 
may be determined from equation (47): 

1 R = 5.4V (47) 

ft a ht 

♦ 

Q R between the gas and surrounding wall of emissivity £ w may 
then be determined from equation (48): 

Q r = £ £ w< r 0 A HT CTg 4 -Tw 4 3 ( 1 hr/3600 sec) (48) 

For gas volumes containing more than one radiating gas, total 
£g is calculated by addition of individual gas emissivities 
and an additional negative correction factor, -A£g, for 
radiant band coincidence which reduces the sum. For internal 
combustion engine exhaust gas products from hydrocarbon 
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